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Abstract 

This  paper  addresses  the  validation  of  an  energy- 
based  inverse  compensator  for  hysteresis  in  ferromag¬ 
netic  transducers.  At  moderate  to  high  drive  levels, 
such  transducers  exhibit  significant  constitutive  non- 
linearities  and  hysteresis  due  to  domain  mechanisms  in¬ 
herent  to  the  constituent  materials.  This  behavior  must 
be  accommodated  to  utilize  the  full  capabilities  of  the 
transducers  in  high  performance  applications.  In  this 
paper,  an  ODE  model  based  on  magnetostatic  energy 
principles  is  employed  to  characterize  the  hysteresis  and 
nonlinear  behavior  of  the  materials.  A  complementary 
ODE  models  the  inverse  behavior  of  the  actuator  and 
is  employed  to  construct  a  filter  which  compensates  for 
the  hysteresis  and  permits  linear  control  design.  The 
capability  of  the  model  and  its  inverse  to  compensate 
for  hysteresis  in  a  magnetostrictive  actuator  is  illus¬ 
trated  through  a  comparison  of  the  input  to  the  inverse 
filter  with  the  actuator  output. 

1.  Introduction 

Ferromagnetic  transducers  provide  the  capability  for 
generating  the  large  forces  and  strains  required  in  a 
variety  of  automotive,  industrial  and  aerospace  appli¬ 
cations.  While  Terfenol-D  is  one  of  the  more  widely 
used  ferromagnetic  core  materials,  the  list  of  materi¬ 
als  currently  employed  for  actuator  design  is  large  and 
growing.  All  of  these  materials  exhibit  a  certain  degree 
of  hysteresis  and  nonlinear  constitutive  behavior  at  the 
high  drive  levels  where  the  materials  are  advantageous 
(see  Figure  1).  In  this  paper,  we  consider  the  experi¬ 
mental  validation  of  a  model  and  its  inverse  which  can 
be  used  to  compensate  for  the  hysteresis  and  nonlinear 
constitutive  behavior. 

To  illustrate  the  goal  in  developing  an  appropriate 
inverse  compensator,  consider  the  prototypical  plant 
depicted  in  Figure  2.  The  input  to  the  actuator  is 
denoted  by  v  while  the  input  which  the  actuator  pro¬ 
vides  to  the  plant  is  denoted  by  u.  Finally,  we  let  Ud 
denote  the  input  specified  to  obtain  the  control  objec¬ 
tive  in  the  absence  of  the  nonlinearities  and  hysteresis 


N(-).  Under  moderate  to  high  drive  conditions,  the 
map  N(-)  contains  sufficient  hysteresis  to  degrade  and 
in  some  cases,  destabilize  linear  control  methods  as  il¬ 
lustrated  in  [8].  To  permit  linear  control  design  through 
inverse  compensation  techniques,  the  desired  input  is 
filtered  through  an  inverse  map  IV_1(-)  or  approximate 
map  IV_1(-)  before  input  to  the  actuator  as  depicted 
in  Figure  2b.  In  this  manner,  one  compensates  for 
the  hysteresis  and  nonlinearities,  and  the  final  signal 
u  =  N(N^1(ud))  is  that  prescribed  by  the  control  law. 

This  technique  for  accommodating  actuator  and  sen¬ 
sor  hysteresis  and  nonlinearities  has  been  employed  in  a 
variety  of  settings,  and  the  reader  is  referred  to  [12,  13] 
for  a  general  development  of  the  theory  in  the  context 
of  adaptive  control  design.  The  models  employed  in  [12] 
to  characterize  the  hysteresis  N(-)  and  inverse  IV_1(-) 
are  based  on  piecewise  linear  Preisach  approximations 
which  are  phenomenological  in  nature.  The  model  em¬ 
ployed  in  this  paper  is  based  on  magnetostatic  energy 
principles  which  facilitates  the  construction  of  the  cor¬ 
responding  filter  and  provides  the  compensator  with  a 
high  degree  of  robustness  with  regard  to  changing  op¬ 
erating  conditions. 

To  motivate  the  theory  used  to  model  the  nonlin¬ 
earities  and  hysteresis,  we  briefly  describe  the  physical 
mechanisms  which  produce  these  phenomena.  Below 
their  Curie  temperatures,  all  ferromagnetic  materials 
exhibit  a  domain  structure  which  consists  of  regions 
where  moments  are  highly  aligned  (domains)  separated 
by  transition  regions  termed  domain  walls.  Hysteresis 
of  the  type  depicted  in  Figure  1  is  typically  attributed 
to  the  restriction  of  moment  rotation  and  domain  wall 
movement  by  inclusions  or  pinning  sites  inherent  to 
the  materials.  This  produces  both  reversible  and  ir¬ 
reversible  components  to  the  magnetization  and  yields 
the  energy  losses  associated  with  Figure  1. 

The  irreversible  and  reversible  domain  processes  are 
modeled  using  the  theory  of  Jiles  and  Atherton  [6,  7]. 
This  quantifies  the  nonlinear  and  hysteretic  relation  be¬ 
tween  the  input  field  and  the  resulting  magnetization. 
The  extension  of  this  theory  to  Terfenol-D  transducers 
was  considered  in  [1,  2,  3,  4]  and  an  inverse  compensator 
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for  Terfenol-D  transducers  was  developed  in  [9].  This 
paper  presents  experimental  validation  of  that  theory. 

The  modeling  techniques  used  to  quantify  the  hys¬ 
teresis  and  nonlinear  constitutive  behavior  are  summa¬ 
rized  in  Section  2  while  the  model  inverse  is  outlined  in 
Section  3.  The  experimental  validation  of  the  inverse 
is  detailed  in  Section  4  for  a  Terfenol-D  magnetostric- 
tive  transducer.  There  it  is  illustrated  that  by  filtering 
a  desired  signal  Ud  by  an  appropriate  inverse  N~l(-) 
before  input  to  the  transducer,  one  obtains  an  input  u 
to  the  plant  which  accurately  approximates  Ud ■  Hence 
the  inverse  can  be  used  to  compensate  for  actuator  non- 
linearities  in  a  manner  which  facilitates  linear  control 
design. 


x  105 


H  (A/m) 


Figure  1.  Relationship  between  the  magnetic  field  H 
and  magnetization  M  for  Terfenol-D  transducer  data. 


(b) 


Figure  2.  (a)  Plant  with  actuator  hysteresis  JV(-); 

(b)  Inverse  compensator  A^_1(-)  for  the  hysteresis. 
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2.  Magnetic  Hysteresis  Model 

The  hysteresis  model  is  based  on  the  theory  of  Jiles 
and  Atherton  [6,  7]  employed  in  the  context  of  magne- 
tostrictive  transducers  [1,  2,  3].  This  theory  is  based  on 
the  tenant  that  ferromagnetic  hysteresis,  of  the  form  de¬ 
picted  in  Figure  1,  is  due  to  the  impediment  of  domain 
wall  movement  by  inclusions  inherent  to  the  materials. 
As  detailed  in  [5],  the  formation  of  domain  walls  at  in¬ 
clusions  or  pinning  sites  is  energetically  favorable  and 
provides  a  minimum  energy  state  for  the  material.  For 
low  field  inputs,  domain  wall  motion  is  reversible  and 
contributes  a  magnetization  Mrev .  At  higher  input  lev¬ 
els,  however,  the  domain  walls  intersect  remote  pinning 
sites  which  corresponds  to  moving  from  one  potential 
energy  well  to  an  adjacent  well.  This  yields  an  irre¬ 
versible  component  to  the  magnetization  Mirr  which  is 
summed  with  the  reversible  component  to  obtain  the 
total  magnetization. 

To  model  the  reversible  and  irreversible  effects,  we 
consider  first  the  hysteresis  free,  or  anhysteretic,  mag¬ 
netization  which  would  result  in  the  absence  of  pin¬ 
ning  sites.  Under  the  assumption  of  constant  stress 
<7o  and  isothermal  operating  conditions,  two  commonly 
employed  models  for  the  anhysteretic  magnetization  are 
the  Langevin  expression 


and  the  Ising  spin  relation 

Man  =  Ms tanh  f— (2) 


(see  [5,  6]).  The  effective  field  given  by 
Hef  f  =  H  +  aM 

incorporates  both  the  applied  field  H  and  interdomain 
coupling  and  magnetoelastic  interactions  through  the 
coefficient  a  and  magnetization  M.  Both  models  are 
derived  by  quantifying  the  probability  of  dipole  mo¬ 
ments  occupying  various  energy  states  through  the  ap¬ 
plication  of  Boltzmann  statistics.  The  higher-order  dif¬ 
ferences  between  the  models  is  due  to  the  underlying 
assumptions  concerning  possible  cell  orientations.  The 
parameters  a  and  a  are  estimated  for  a  given  transducer 
through  a  least  squares  fit  to  the  data. 

As  detailed  in  [6],  the  irreversible  magnetization  is 
quantified  by  determining  the  energy  required  to  re¬ 
orient  dipoles  in  the  presence  of  an  applied  field.  This 
yields  the  differential  equation 


dMirr  _ ^  Man  Mir  r 

dH  kS  —  a  ( Man  —  Mirr ) 


(3) 


where  the  parameter  k  provides  an  average  measure 
of  the  energy  required  to  re-orient  domains  and  hence 


translate  domain  walls.  The  parameter  8  is  defined  to 
have  the  value  +1  when  dH  >  0  and  —1  when  dH  <  0 
to  guarantee  that  pinning  always  opposes  changes  in 
magnetization.  In  applications,  5  can  be  directly  de¬ 
termined  from  the  magnetic  field  data.  Finally,  the 
parameter  5  defined  by 


f  0 


5={ 


{dH  >  0  and  M  >  Man} 
or  {dH  <  0  and  M  <  Man} 


(,  1  ,  otherwise. 


incorporates  the  physical  observation  that  following 
field  reversal  at  the  loop  tip,  the  changes  in  magne¬ 
tization  are  reversible  until  the  anhysteretic  value  is 
reached  [7]. 

The  second  component  of  the  magnetization  is  the 
reversible  magnetization  which  reflects  the  degree  to 
which  domain  walls  bend  before  attaining  the  energy 
necessary  for  translation.  As  derived  in  [6],  to  a  first 
approximation,  the  reversible  magnetization  is  given  by 

Mrev  =  c(Man  -  Mirr) .  (4) 


The  reversibility  coefficient  c  can  be  estimated  from  the 
ratio  of  the  initial  and  anhysteretic  differential  suscep¬ 
tibilities  [7]  or  through  a  least  squares  fit  to  data. 

The  total  magnetization  is  then  given  by 


A4.  —  Mrev  +  Aljrr 


(5) 


with  Mirr  and  Mrev  defined  by  (3)  and  (4)  and  the  an¬ 
hysteretic  magnetization  given  by  (1)  or  (2).  To  provide 
a  relation  which  facilitates  inversion,  it  is  advantageous 
to  express  the  output  magnetization  as  a  function  of  the 
input  field.  When  the  Ising  spin  model  is  used  to  char¬ 
acterize  the  anhysteretic  magnetization,  the  expressions 
(2),  (3)  and  (4)  can  be  consolidated  to  yield 


dM 

~dH 


.  .  dMirr  dMan 

(1  -c)  — rrr-  +c 


5(1  -c) 


dH  dH 

Man  ~  M^r 


+  C- 


dMn 


kd  -  a(Man  -  M^r)  '  dH 
8  [Mstanh  ( H+“M )  -  M] 


kd  -  a  [Ms tanh  ( H+(lM )  -  M] 


cMs  2 
H - sech 


H  +  aM 


1+0; 


dM 

~dH 


where 

T{H ,  M)  = 


1  - 


cMsa 


sech2 


8  [Ma tanh  -  M] 

k8  -  a  [Ma tanh 


cMs  ~(H+aM 
H - sech 


If  one  employs  the  Langevin  expression  (1)  rather  than 
the  Ising  spin  relation  (2)  for  the  anhysteretic,  the  func¬ 
tion  T  is  given  by 


T{H,  M)  = 


1 


1+CMso|#(M) 

8[MsZ{S+2M.{  -  m] 

k8  -  a  [MSC  -  M] 

dC  (H  +  aM' 


(8) 


-cM, 


dH 


where  the  Langevin  function  is  defined  by 
C(z)  =  coth(z)  -  -  . 


The  relation  (6)  with  T  given  by  (7)  or  (8)  is  then  em¬ 
ployed  when  quantifying  the  actuator  hysteresis  N(-). 
In  both  cases,  the  parameters  8  and  8  can  be  computed 
directly  from  the  measured  field  data  or  the  computed 
values  of  the  magnetization  and  anhysteretic  magneti¬ 
zation.  The  saturation  magnetization  Ms  can  often  be 
determined  directly  from  the  measured  data  or  a  pri¬ 
ori  knowledge  of  the  material  behavior.  Hence  only  the 
parameters  a,  a,  c  and  k  must  be  determined  through  a 
least  squares  fit  to  data  or  adaptive  parameter  estima¬ 
tion  techniques. 

Finally,  we  note  that  the  time-dependent  dynamics 
of  the  magnetization  can  be  specified  through  the  chain 
rule 


dM 

dt 


T{H,  M) 


dH 

dt 


This  relation  should  be  employed  only  for  quasistatic 
processes,  however,  since  the  model  (6)  neglects  the 
elastic  properties  of  the  material  as  well  as  time- 
dependent  loss  mechanisms. 


where  a  = 
level  is  then 
equation 


The  magnetization  at  a  given  field 
specified  by  the  solution  to  the  differential 


dM 

~dH 


F(H,  M) 


(6) 


M(H0)  =  M0 


3.  Inverse  Hysteresis  Model 

The  hysteresis  model  in  the  form  (6)  is  amenable  to 
inversion  through  consideration  of  a  complementary  dif¬ 
ferential  equation.  For  actuators  in  which  the  anhys¬ 
teretic  is  modeled  by  the  Ising  spin  model,  the  inverse 
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is  specified  through  solution  of  the  differential  equation 

dM -1  _  1 

~dir  ~  f(m-\h)  (9) 

M-\H0)  = 

with  T  defined  by  (7).  Similarly,  employment  of  (8) 
provides  the  inverse  when  the  Langevin  model  is  used 
to  characterize  the  anhysteretic  magnetization. 

In  the  control  nomenclature  of  Section  1,  the  defini¬ 
tion  (9)  can  be  combined  with  passive  actuator  or  struc¬ 
tural  components  to  provide  the  exact  inverse  77_1(-)  if 
the  parameters  Ms,  a,  a,  c  and  k  are  known  exactly.  For 
the  examples  in  this  paper,  it  is  assumed  that  suitable 
parameters  have  been  obtained  through  a  least  squares 
fit  to  data  and  the  exact  inverse  is  employed.  Simi¬ 
larly  it  provides  an  approximate  inverse  IV_1(-)  which 
is  suitable  for  adaptive  parameter  estimation  and  con¬ 
trol  methods  if  the  parameters  are  unknown  or  slowly 
change  due  to  changing  operating  conditions.  Details 
regarding  the  design  of  control  methods  which  employ 
an  approximate  inverse  can  be  found  in  [12]. 

4.  Experimental  Validation 

To  illustrate  the  performance  of  the  model  and  in¬ 
version  technique,  we  consider  the  characterization 
and  hysteresis  compensation  for  a  prototypical  mag- 
netostrictive  transducer  employed  at  Etrema  Products. 
As  depicted  in  Figure  3,  the  primary  components  in  the 
transducer  are  the  driving  Terfenol-D  rod,  a  wound  wire 
solenoid  which  surrounds  the  rod,  and  a  prestress  mech¬ 
anism  comprised  of  Belleville  washers  which  maintain 
the  rod  in  a  state  of  compression.  The  results  reported 
here  were  obtained  with  a  prestress  of  approximately 
Go  =  1.0  ksi.  The  Terfenol-D  (Tbo.3Dyo.7Fe1. 92)  rod 
had  a  length  of  0.254  m,  a  diameter  of  0.020  m,  and 
was  composed  of  four  laminae  to  minimize  eddy  cur¬ 
rent  losses.  The  drive  solenoid  was  constructed  from 
14  AWG  round  wire  with  an  air  gap  left  between  the 
coil  and  Terfenol  rod  to  permit  cooling.  Additionally, 
a  pickup  coil  comprised  of  34  AWG  wire  was  wrapped 
directly  around  the  rod  for  sensing  purposes. 

For  characterization  and  validation  of  the  inverse 
compensator,  it  was  necessary  to  measure  the  field  H 
and  magnetization  M  produced  in  the  transducer.  The 
current  I(t )  was  used  to  compute  the  field  using  the 
relation  H(t )  =  nl(t)  where  n  is  the  number  of  turns 
in  the  drive  coil  per  unit  length.  We  note  that  this 
relation  ignores  end  effects  in  the  rod  and  the  reader 
is  referred  to  [2,  3]  for  details  concerning  the  construc¬ 
tion  of  an  empirical  filter  which  incorporates  the  spatial 
behavior  of  the  rod.  From  the  induced  voltage  in  the 
pickup  coil,  the  Faraday-Lenz  law  was  used  to  compute 
the  temporal  derivative  of  the  magnetic  induction  B. 
Integration  then  yielded  the  induction  and  magnetiza¬ 


tion  M  =  —  H.  Finally,  we  note  that  the  data  was 

collected  at  1  Hz. 

The  first  step  in  the  construction  of  the  inverse  com¬ 
pensator  was  the  determination  of  parameters  a,  a,  k,  c 
and  Ms  so  that  the  model  accurately  characterizes  the 
transducer.  While  initial  values  for  the  parameters  can 
be  directly  determined  from  measured  values  of  the  dif¬ 
ferential  susceptibility  and  field  [7],  or  a  priori  manu¬ 
facturer  specifications,  final  values  are  typically  deter¬ 
mined  through  a  least  squares  fit  to  the  data. 

The  model  response  obtained  with  the  parameters 
Ms  =  7.65  x  105  A/m,  a  =  1.5  x  105  A/m,  a  =  0.032, 
k  =  5000  A/m,  and  c  =  0.18  is  compared  with  the 
data  in  Figure  4.  We  note  that  while  the  data  exhibits 
changes  of  concavity  in  the  burst  region  which  are  not 
incorporated  in  the  model,  the  model  accurately  quanti¬ 
fies  the  hysteresis  and  saturation  behavior  of  the  trans¬ 
ducer. 
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Figure  3.  Cross  section  of  the  Terfenol-D  magne- 
tostrictive  transducer. 
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Figure  4.  Experimental  M-H  data  ( - )  and  model 
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The  inverse  (9)  with  the  above  parameter  values  was 
then  employed  as  a  filter  for  the  physical  transducer  in 
the  manner  illustrated  in  Figure  2.  A  1  Hz  triangular 
input  Ud  was  applied  to  the  inverse  to  generate  a  signal 
v  which  was  input  to  the  transducer.  The  output  u 
from  the  transducer  was  then  measured  to  determine 
the  accuracy  of  the  compensation  technique. 

The  input  signal  Ud  is  compared  with  the  transducer 
response  u  in  Figure  5.  It  is  observed  that  while  u 
exhibits  slight  deviations  due  to  unmodeled  domain  ro¬ 
tation  effects  in  the  burst  region,  the  overall  fit  is  ac¬ 
curate.  Hence  the  inverse  provided  by  this  model  ade¬ 
quately  compensates  for  the  hysteresis  and  saturation 
behavior  of  the  transducer. 

From  the  perspective  of  control  design,  it  is  necessary 
to  quantify  the  degree  to  which  the  transducer  can  ac¬ 
commodate  various  inputs.  Aspects  of  this  question 
are  illustrated  in  Figure  6.  Here  the  field  v  prescribed 
by  the  inverse  compensator  is  compared  with  the  field 
which  is  actually  measured  in  the  transducer.  The  two 
agree  through  the  low  and  middle  drive  ranges  but  the 
measured  field  does  not  attain  the  peak  values  due  to 
limitations  in  the  hardware  and  amplifier.  This  can  be 
alleviated  by  limiting  the  rate  and  magnitude  of  pre¬ 
scribed  inputs. 

5.  Concluding  Remarks 

The  results  presented  here  illustrate  both  the  perfor¬ 
mance  of  a  domain  wall  model  for  quantifying  hysteresis 
and  nonlinear  constitutive  behavior  in  a  magnetostric- 
tive  transducer  and  the  development  of  an  inverse  com¬ 
pensator  based  on  this  model.  The  model  and  its  in¬ 
verse  are  formulated  as  ODE  which  are  implemented  us¬ 
ing  low-order  approximation  techniques.  The  efficiency 
of  the  method  is  further  augmented  by  the  fact  that  the 
model  and  inverse  rely  on  only  five  parameters,  some 
of  which  can  be  directly  approximated  from  a  priori 
specifications  or  measured  properties  of  the  data. 

A  comparison  of  the  input  Ud  to  the  model  inverse 
with  the  output  data  u  from  the  physical  transducer 
illustrates  that  the  technique  accurately  compensates 
for  the  hysteresis  and  nonlinear  behavior  of  the  trans¬ 
ducer.  Furthermore,  the  observation  that  unmodeled 
phenomena  are  not  exacerbated  by  the  inverse  hints  at 
the  stability  of  the  technique. 

The  design  of  linear  control  methods  which  employ 
inverse  compensators  has  been  well  established  in  the 
literature  (e.g.,  see  [12])  and  initial  numerical  investi¬ 
gations  have  indicated  the  viability  of  LQR  designs,  for 
structures  employing  magnetostrictive  transducers,  us¬ 
ing  this  compensation  technique  [9] .  Current  investiga¬ 
tions  are  focussed  on  the  experimental  implementation 
of  linear  control  methods  which  utilize  this  compen¬ 
sator,  and  the  extension  of  these  techniques  to  ferro¬ 
electric  materials  [10,  11]. 


Figure  5.  The  input  Ud  to  the  inverse  ( - )  and 

output  u  from  the  transducer  ( - ) . 


Figure  6.  The  field  v  prescribed  by  the  inverse  as 

input  to  the  actuator  ( - )  and  the  field  measured  in 

the  actuator  ( - ). 
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